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Numerical evidence for a finite-time singularity in ideal 
3D magnetohydrodynamics (MHD) is presented. The simula- 
tions start from two interlocking magnetic flux rings with no 
initial velocity. Curvature shrinks the rings until they touch 
and current sheets form between them. The evidence for a 
singularity in a finite time t c is that the peak current density 
behaves like ||J||oo ~ l/(t c — t) for a range of sound speeds 
and initial conditions. For the incompressible calculations 
||<^||oo/||J||oo — * const. In resistive reconnection the magnetic 
helicity is nearly conserved while energy is dissipated. 



To date it is not known whether or not the ideal mag- 
netohydrodynamic (MHD) and Euler equations are reg- 
ular. Regularity means that for arbitrary smooth initial 
data the velocity u and/or magnetic field B and all of 
its derivatives remain finite for all times. If regularity is 
not established, one cannot rule out the possibility of a 
finite time singularity, which for ideal MHD implies sat- 
isfying a new mathematical constraint that the vorticity 
u> = V x u and current density J = V x B must obey 



occur in a finite time t — t c in any quantity (e.g. in any 
derivative of u, however high) unless 



I(t c )= / [HojIUC*) + ||J||oo(*)] rf* ^ oo, 



(1) 



for there to be a singularity. Singularities associated with 
shocks are also possible, but are not considered here, be- 
cause they would not occur in the incompressible limit. 
The physical significance of a finite-time singularity is 
that fast vortex and/or magnetic field reconnection could 
be possible once viscous and resistive effects are restored. 
A theoretical demonstration of a mechanism for fast mag- 
netic reconnection would be significant in a variety of 
problems in plasma physics including experimental stud- 
ies related to magnetically confined fusion 0, dynamos, 
the earth's magnetic field and the solar corona, where 
reconnection is associated with flares HQ and coronal 
heating ||||. This letter will show preliminary evidence 
for a possible singularity for ideal MHD using numeri- 
cal simulations of three-dimensional linked magnetic flux 
rings and emphasize its relationship to fast magnetic re- 
connection. 

Significant progress has been made recently for the Eu- 
ler case where a numerical study of a pair of antiparallel 
vortex tubes has produced strong evidence for the for- 
mation of a singularity in the Euler equations in a finite 
time . The numerical evidence was consistent with an 
analytic constraint for Euler 8 that no singularity can 



Htc) 



(2) 



Here ||...||oo is the L°° norm, or maximum, in space. This 
theorem shows that when searching for a singularity in 
the 3D Euler equations, the only quantity that needs to 
be monitored numerically is Hu^loo. Furthermore, if a 
singularity of the form ||w||oo ~ (t c — i)~ 7 is observed in 
a numerical experiment then 7 must obey 7 > 1 for the 
observed singular behavior to be genuine. The numerical 
comparisons M found 7 = 1, a scaling that is consistent 
dimensionally in addition to being consistent with (|^). 
The generalization to MHD is @). 

An important feature of the analysis of the Euler 
calculations that should be used when analysing other 
flows with possible singularities is that singular behav- 
ior should be demonstrated by several independent tests. 
One set of tests predicted mathematically is that in addi- 
tion to the l/(t c — t) 1 behavior of ||w||oo> all components 
of || Vu||oo should diverge as l/(t c — t) 1 Again, 7 = 1 
is expected. Another test that was found numerically |?]] 
is that the rate of enstrophy production is 



dn/dt 



uJiCijUj dV ~ l/(t c - t) 



(3) 



where enstrophy is ^fi = J \u)\ 2 dV. 

The question of fast magnetic reconnection has been 
addressed using both time evolving calculations and 
steady-state analyses. In two-dimensional calculations, 
fast reconnection is inhibited by material that cannot 
escape easily from between thin current sheets. This 
led to early suggestions that a singularity in 2D MHD 
is precluded by the formation of current sheets [jlO| and 
is consistent with a recent result that there can be no 
singularity at 2D nulls JTl[| . New experimental evidence 
designed to test two-dimensional steady-state theo- 
ries shows a slow resistive timescale Jl2|-p^| instead of a 
faster timescale resulting from slow shocks |l5| l and re- 
lated theories |l6)| . In three dimensions, the extra degree 
of freedom could allow material to escape more readily 
from the current sheets. This has prompted reconnec- 
tion simulations in three dimensions, for example starting 
with orthogonal, straight flux tubes in pressure equilib- 
rium and a small velocity field to push them together 
]l7t . Explosive growth of ||J|| has been seen before [ [l8| , 
but not in the context of singularities. 

We consider the equations for an isothermal, compress- 
ible gas for a given sound speed c s with a magnetic field 
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of the form 
du 



dt 



u-Vu-c 2 Vlnp+ J X B +- (V 2 u+ ±VV • u) , 

p p K J ' 



91np 
dt 



dA 
~dt 



u ■ V In p — V ■ u, 



u x B + ?7V 2 A, 



(4) 



where B = V x A is the magnetic field in terms of the 
magnetic vector potential A, u is the velocity, and p 
is the density. In the ideal limit, the resistivity r\ and 
the viscosity p are set to zero. The magnetic field is 
measured in units where the permeability is unity. Pe- 
riodic boundary conditions are adopted in a domain of 
size (27r) 3 . Our time unit is the sound travel time over a 
unit distance. The equations are advanced in time using 
a variable third-order Runge-Kutta timestep and sixth 
order explicit centered finite differences in space. 

The equations for the incompressible case are the 
same except that the equation for p is replaced by the 
divergence-free condition on velocity V ■ u to determine 
the pressure. The magnetic field B, rather than the vec- 
tor potential A, is used as a fundamental variable. The 
equations are advanced in time using a spectral coloca- 
tion method with the 2/3-rule and variable third-order 
Runge-Kutta timestep. The maximum resolution used 
384 3 mesh points. 

The incompressible, ideal MHD equations conserve 
total energy E = | /(u 2 + T5 2 )dV, magnetic helicity 
Hb = J A • B dV and the cross helicity He = / u • B dV. 
The helicities can be used to describe aspects of the topol- 
ogy (l^] . Since Hg has one less spatial derivative than the 
energy, spectrally it should dissipate more slowly than 
energy |20| when p, v ^ 0. 

The initial conditions used in the present study will all 
be of two linked, magnetic flux rings. This condition is 
chosen because it has the advantage that no velocity field 
needs to be imposed in the initial conditions since the 
tension from the curvature of the rings induces a velocity 
by shrinking the rings. Due to the initial linkage, it also 
yields nearly maximal Hb and therefore is an excellent 
choice for studying the dissipation of Hb versus energy 

A variety of different ring thicknesses and angles be- 
tween the initial rings have been investigated. This letter 
will discuss in detail only cases where the rings are or- 
thogonal and are just touching. There are three distances 
that determine the initial condition that evolved into the 
structures in Figure |l|: the radii of the rings R, the thick- 
ness r where the flux goes smoothly to zero, and the sep- 
aration of their centers from the origin A. The separation 
of the rings is 2(R— A). An initial profile across the ring 
that gives \B\ = 1 in the center and goes smoothly to 
\B\ = at r = r a is taken from the Euler calculations 
0. R = 1 and A = 0.5 for all the cases, so that they go 



through each other's centers. The initial condition for the 
compressible calculations to be reported used r = 0.5. 
r = 0.5, 0.65 and 0.8 for the ideal incompressible cases. 
Following the example from the Euler case 0, the fol- 
lowing hyperviscous filter was applied in Fourier space to 
the initial condition only: exp[— (k/k max ) ), where & max 
is between 14 and 20. As a result, the maximum initial 
magnetic field Bq is slightly less than one. For the com- 
pressible calculations, c s is varied between 0.1 and 10, so 
the initial plasma beta, /3q — 2c 2 /i?g, varied between 0.02 
and 200. Only nearly incompressible flo = 2 calculations 
are presented. The initial density was uniform and unity. 
While what is being presented is nearly incompressible, 
preliminary low plasma beta calculations as well as com- 
pressible vortex reconnection pl| suggest that compress- 
ibility might enhance reconnection rates and singularity 
formation. 

All of the compressible, resistive calculations were 
straightforward runs from t = with viscosities and re- 
sistivities chosen for a given resolution. The strategy 
for reaching the highest possible resolution for the ideal 
calculations followed the example of how a possible sin- 
gularity in incompressible Euler was demonstrated [Q. 
First, the ideal calculations do not contain any numeri- 
cal smoothing and have been run only so long as numeri- 
cal dissipation was insignificant. This approach was used 
because experience has shown that artificial smoothing 
results in artificial dissipation which can obscure the dy- 
namics of the ideal case. The lower resolution calcula- 
tions were started from t = 0. Since it would be too ex- 
pensive to run higher resolution calculations from t = 0, 
they are being remeshed at intermediate times. 

There was no initial velocity field. The first, and short- 
est, phase after initialization was that due to the curva- 
ture of the flux rings their diameter R shrinks. This 
automatically brings the flux tube rings into contact and 
the current of one ring begins to overlap the magnetic 
field of the other. This is necessary for the Lorentz force 
to be significant and for a strong interaction to begin. 

To give an overall view of the flow, Figure |l| shows 
the three-dimensional structure from a j3 = 0.5, 192 3 
resistive, compressible calculation, just before and after 
the estimated singular time, t = 2 shows nearly ideal 
evolution from the initial motionless, perfectly tubular 
flux tubes. The dominant feature is the indentation in 
each ring. This is because in the cores of the flux rings 
the field is strongest and the magnetic curvature force 
largest. The region going singular appears in Figure 2a. 
as twisted current sheets within the center of Figure 1 a. 
A surprising property of the reconnection process is that 
slices show that the size of the entire structure shrinks. 
Some perspective on this can be obtained by comparing 
the structures at t = 2 and t = 3. 
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FIG. 1. Resistive calculation using compressible code. 
t — 2 shows evolution during the nearly ideal phase, t — 3 
shows a partially reconnected state with Hb converted more 
into new twist between the remnants of the original tubes 
than into writhe within reconnected tubes. 



By t = 3 some reconnection has occurred. Magnetic 
helicity Hb is nearly conserved, decreasing linearly at a 
very slow pace as it is converted into writhe |2^| or new 
twist, and energy is dissipated more rapidly. A detailed 
examination of this structure will be studied elsewhere. 




FIG. 2. Three-dimensional level surfaces of the magnitude 
of the current density in the resistive calculation, t = 2 shows 
a double saddle surface, which at t — 3 is broken up into two 
disjoint pieces. 



In a resistive calculation, the dissipation is concen- 
trated on a current sheet that forms where the two flux 
rings come into contact. Level surfaces of the current 
density near its peak value (Fig. [2]) show a twisted, sad- 
dle shaped, double-sheet structure before the estimated 
singular time, which separates into two disjoint sheets 
centered around the points of maximum current density 
after the singular time. 



Figure I shows 1/||J|U and l/(||J||oo + IMU) for 
the resistive calculations. There is a strong tendency 
in favor of linear behavior similar to that observed for 
3D Euler 0. Extrapolating from before t — 2.3 to 
1/Jmax — suggests that t c ns 2.5. For c s > 0.5, that 
is more incompressible, roughly the same singular time 
would be predicted. For c s < .5, that is more compress- 
ible, different behavior is indicated, but the trend toward 
1/ 1| J || 00 ~ (t c — t) remains. 
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FIG. 3. Evolution of l/||J||oo and 1/(||J|U + ||w||co). 
r — 0.5 for /3q = 2. Solid and dashed lines refer to ideal cal- 
culations with different resolution and the dotted line is for 
a resistive calculation with r) = fJ, = 0.002. The dash-dotted 
line gives a linear fit to the data. 

For Euler, recall that stronger evidence for a possible 
singularity was obtained by monitoring the L°° norm for 
an additional strain term and the enstrophy production, 
as well as |[u>[[oc @- Following that reasoning, we need to 
know the behavior of two first derivatives of the magnetic 
or velocity fields plus a global production term. There- 
fore, we propose looking at the behavior of ||J||oo, ||<*>||ooj 
and the production of Q u+ j = \ J (|u>| 2 + | J| 2 ) dV, 

Puj = /(W - LOidijJj + 2e ijk J idjt e ik ) dV, (5) 

where ey = h(ui,j + and = \{Bij + Bj^) are 
the hydrodynamic and magnetic strains. The terms in 
Pnj, in order, are the vortex stretching already known for 
Euler, a new vorticity production term and a new current 
producing term. All three tests should go as l/(t c — t) 
once sufficiently singular solutions are obtained, which we 
have not yet achieved. Therefore, the present objective 
is trends in the right direction, rather than conclusive 
evidence for the existence of a singularity. 

Figure || plots the three proposed tests for the r = 
0.65 incompressible case. The 384 3 run was used for the 
l/||J||oc and l/H^lloo tests for t > 1 and the 192 3 run for 
1/Pqj for all time. 
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FIG. 4. Comparison of the evolution of 1/ 1| J |] oo , 1/||w||do 
and Pnj for r a = 0.65. The inset shows || oo / 1| J || oo 
for r =0.5, 0.65 and 0.8, all of which are converging to 
1 1 <-«-' 1 1 oo / 1 1 J||oo ~ 0.5, demonstrating possible self-similarity. 

In the case of Euler it was shown that ratios of all 
L°° norms of first derivatives were approaching constant 
values Q, which would be consistent with self-similar 
behavior near the point going singular. Therefore we ex- 
pect that || w || do /|| J||oo should approach a constant value 
here. Initially the velocity, and vorticity, are zero, so 
these must build up before behavior for ||w||oo and P07 
consistent with l/(t c — t) appears. The inset in Figure [| 
shows that all three cases appear to be converging to a 
value of ||w||oo/||J||oo ~ 0.5. r = 0.65 converges the 
soonest, near t = 1.6 and together with the convergence 
of ||w||ao/||J||oQ — *■ 0.5 is our best evidence to date for 
a possible singularity in ideal MHD. In addition, 1/Pqj 
appears to join the same (t c — t) near t = 1.9. 

While 1/ 1| J || 00 for the r = 0.5 incompressible case 
is comparable to the compressible case and shows the 
largest range of 1/ 1| J || 00 ~ {tc — t) behavior, it is not 
strong evidence for a singularity because the inset in Fig- 
ure f| shows that its l/Hu^oo up to t = 2.0 does not have 
the same value for t c . However, note that the production 
term for J 2 in (||), 2£ijkJidj^eik, does not involve vortic- 
ity, but only strain terms. Unlike Euler, vorticity seems 
to be playing only a secondary role for ideal MHD. It is 
not essential for vorticity to blow up for the compressible 
cases, and it does not. Related to this, the current sheets 
near ||J||oo m Figure || are twisted. This also suggests 
that similar hydrodynamic initial conditions should be 
revisited to determine their reconnection rates p3| . Fur- 
ther analysis of local production terms in proposed 600 3 
calculations should help answer these questions. 

In conclusion, we have presented numerical evidence 
for a finite-time blow-up of the current density in ideal 



MHD in the case of interlocked magnetic flux rings. In 
the resistive case one would not expect there to be a 
singularity. Instead, arbitrarily thin current sheets will 
form, depending on how small the resistivity is 0] . This 
can lead to significant dissipation whose strength is vir- 
tually independent of resistivity. The astrophysical sig- 
nificance of current sheets as a consequence of tangential 
discontinuities of the field has also been stressed in a 
recent book by Parker However, the possibility of 
a finite time singularity in ideal MHD is not commonly 
discussed in connection with fast reconnection. 
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